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Abstract 



A problem of monoenergetic particles pulse reflection from half-infinite 
stratified medium is considered in conditions of elastic scattering with ab- 
sorbtion account. The theory is based on multiple scattering series solu- 
tion of Kolmogorov equation for one-particle distribution function. The 
analytical representation for first two terms are given in compact form 
for a point impulse source and cylindric symmetrical detector. This pub- 
lication is an authorized translation (by S. Leble) of a bad-published paper 



from 20.06.1980 (available in its original version: http://www.mif.pg.gda.pl/krrizm/page/leble/publications.html 

140. A. Buzdin, S. Leble Lidar Problem Solution in Double-Scattering 

Approximation VINITI N2536-80dep (1980) 126 [PDF]). Reading recent 

articles on the LIDAR sounding of environment (e.g. Atmospheric and 

Oceanic Optics (2010) 23: 389-395, Kaul, B. V.; Samokhvalov, I. V.) one 

recovers standing interest to the related direct and inverse problems. A 

development of the result for the case of n-fold scattering and polarization 

account as well as correspondent convergence series problem solution of 

the Kolmogorov equation will be published in nearest future. 

1 Introduction 



A problem of monoenergetic particles pulse reflection from half-infinite strat- 
ified medium is considered. The particles scattering is supposed elastic, its 
interaction is negligible. Such problem relates a LIDAR probing of atmosphere. 
It is supposed that such process is described by Kolmogorov equation [1 for 
one-particle distribution function. 

The expansion of the distribution function vifith respect to scattering mul- 
tiplicity allows to obtain a closed expression for the series terms as multiple 
integrals. 
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The main purpose of this paper is a simplification of analytic expressions 
and investigation of its applicability range in terms of geometric parameters and 
macroscopic cross-sections. The point impulse source and cylindric symmetrical 
detector is considered. 

Computations by Monte-Carlo method show a quick convergence of the scat- 
tering multiplicity expansion hence the one- and double scattering contributions 
are enough for the practical LIDAR response evaluation. 

2 The problem formulation 

The equation for the probability density / = f{t,f,v,t') has the form: 

^[dt'f + v-Vf] = -at{z)f- J aicos-f,z)fdn, (1) 

Where t'-timc, dO - solid angle; cr— bulk differential cross-section of elastic 
scattering to the angle j; at - the sum of / adSl and total cross-section of 
absorbtion, and 

V 

- = {sin9cos<p,sin6sin^,cos6). 

In spherical coordinates: r, 6, (j) the scattering angle is expressed as 

cos'y = cos6cos6' + sin0sin6'cos{(f) — (j)'). (2) 

We suppose that the scattering is elastic, \v\ does not change while the scattering 
process occur. Initial conditions are represented by distributions 

f(Q,x,y,z) = V5{x)5{y)5{z)5{6). (3) 

It means that we built a solution for the probability density as a weak limit 
(when t ^Qi) to 5 — function at f > 0. The distribution 5{6) is chosen as 

{5{e),^{e,<t>)) = r i;{0,<p)d<p. (4) 



3 Solution 

Let us denote t = ct'. A solution is searched as an N-fold scattering expansion 

/ = /o + /l + /2 + ... (5) 

We choose for /o, 

Lfo = -K— + sinOcosd)^^ + sinOsincb^^ + cosO'^^ = —crt{z)fo, (6) 
at ox ay oz 

and initial condition, 

/o(0, X, y, z) = V5{x)S{y)S{z)6{0). (7) 
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For n > 1, the expansion coefficients are defined by 

Lfn+i = -ai{z)fn+i + / / a2{cosj,z)fnsin0'd6'd(l)' . (8) 
Jo Jo 

with initial conditions for n > 

/n+l|t=0 = 0. (9) 

Such expansion is usually named as expansion by folds scattering. The expres- 
sion for /o is easily found: 

/o = ^S{x)S{y)S{z - t)6{e) exp[- /* at{z - T)dT]. 

^TI" Jo 

Let us denote a function E via 

E{t, z, 6) = exp[- / at{z - TC0s6)dT], (10) 
Jo 

then for fn+i we obtain the reccurence 

fn+i{x,y,z,e,(f),t) = 
Jo E{t, z, 9) J^"" a{cos-f, z - tcos9) (11) 
fn{x — TCOsOcoscp, y — Tsin9sin(j), z — tcos9, 9', 0', t — T)sin9' d9' d(j)' dr. 

Particulary, for the n = 1, it yields 

fi= E{t, z, 9)E{t — T,z — TC0s9, 0)a{cos9, z — tcos9) 
Jo 

S{x — Tsin9cos^)S{y — Tsin9sin^)6{z — tcos9 — {t — T))dT 
and for n — 2 

f2 = /o E{t2, Z, 9) /q^" /q" C7(cOS7, Z - T2C0S9) /q E(ti,Z - T2C0S9, 9') 
E{t — T2 — Ti, Z — T2C0S9 — T\COs6' , 0)a{cos6' , Z — T2C0S6 — T\COs9') 

5{x — T2sin9cos<j} — TisinO' cos(j)')5{y — T2sin9sin(j) — Tisin9' sincj)') 
5{z — T2Cos9 — Ticos9' — {t — T2 — Ti))dTisin9' d9' d<p' dT2. 

(12) 

It is convenient to rcrpresent the integrands as distributions in x, y, z space 
parametrized by 9, 0, ti,T2,9i,4>i. Integrations by n, r2, 0i, (/>i is understood as 
integration of the distribution by these parameters. For example, /i acts as on 
a function tp from Shwartz space as 

{fi{t,x,y,z,9,(t)),'il){x,y,z)) = 
E{t, tcos9 + t-T, 9)E{t -T,t-T, 0)a{cos9, t - t) 
ip{Tsin9cos(j), Tsin9sin(j), tcos9 + t — T)dT. 
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4 Number of particles rate 

Our next problem is evaluation of number of particles which enter the round 
area of radius po la-ying in the plane xy with center in the origin and having 
velocity vectors inclined to z-axis within the angle interval 6 ^ [t: — Oq^-k]. The 
angle relates the aperture angle of a receiver ( e.g. lidar telescope window). 
By its direct sense the number is proportional to 

I{t) ^lim— / {f{x,y,z,e,(j),t),ip{x,y,z))sined(j)de, (13) 

^* JjT-eo Jo 

where ip{x,y,z) = 1 for internal points of the domain + j/^ < q < 

z < At\cos6\ and zero outside. 

Contribution of n — fold scattering /„(i) is obtained by substitution /„ 
instead of / in the Eq. p^ . Let us evaluate one- and two- fold scattering. 

liniAt ^ £_g^ Jo'' E{t, Tcose + t-T, e)E{t -T,t-T, O)cr{cos0, t - r) (14) 
il)(Tsin9, 0, TcosO + t — T)dTsin6 



The independence of the integrand on (/) is taken into account. 

In the problems of LIDAR probing the aperture angle is small 0o << 1, 
therefore we estimate cos0 as —1. Introduce new variables z — t — t,^ = 
TsinO. Intersection of the integration domain in (jl4p with the area in which -0 
is nonzero, satisfy the following inequalities (in new variables) 

< 2z - i < At, (15) 

0<e<e(t-^), 

where e = tgO^. When 5 > the third inequality from ((T5|) follows from the 
first ones, hence 



h{t) = 

limAt^o C E{t - z, 2z - t, 7t)E{z, z, 0)a(~l, z)^ 



(16) 



'^(;^,0,^)i?(-,-,0)<7(-l,-), (17) 



^2 ' y 2' ' ' '2' 

which results in LIDAR formula 



'°eM~ at{z)dz]o{-l,-). (18) 







t2 

Contribution of the two-fold scattering ([T2|) gives 



limAt^o^/^_e^/J'/o/o ^/o''Ecr(cos7,Ti (0056*1 - 1) + i - T2)cr(cos0i, +t - n - T2) 
'ip{Tisin9i + T2sin92C0S(j), T2sin92sin(j), ticos9i + T2Cos92 +t — ti — T2)d9idT2dTid(f)d92 

(19) 
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where 

E = exp[- Jp ""^ at{t -T2-T1- T[)]dT[ JJ"' a-t((ri - r')cos6'i + i ~ T2 - Ti)dT' - 

Jq^ <7t{TlCOsdl + {t2 — t')cOs62 + < — T2 — Tifdr' 

(20) 

It was used that the integrand depends only on (j) — (j)i~(f>2 and, due to definition 
of ■0 and C0S7, is even in (j). 

Let us introduce new variables: 

Zl = i - T2 - Ti , 
Z2 = TCOsOl +t — Tl — T2, 

^=Tisinei, (21) 
= T2sin92, 

P = [t^ sin6'^ + 2TiT2sin9isin92 + r'^sinO'lY/'^ 

The sense of 2:1.2, P has transparent intuitive meaning: zi.2 are heights of scat- 
tering points, p is the polar coordinate of a trajectory end. In the domain of 
consideration there is a one-to-one correspondence between old and new coor- 
dinates. The domain is fixed by ^ > 0, 77 > 0, and 



n = \J{zi - Z2Y T2^t - Zi- Tl, 



cosOi 



^{zi-z^r+i-' (22) 



The Jacobian evaluation gives: 

Jsin9isin92 = 7 ^ „ , . , , (23) 

(T1T2) \cos9i\sin(p 

that is obtained via transformation formulas. Below, as in derivation of the 
one- fold scattering, we shall approximate cos9i as 1. 

Let us obtain conditions that restrict the integration domain taking into 
account the definition og ip. The conditions for p are obvious. From the problem 
statement it follows also that zi > 0, Z2 > 0. Due to |cos(/)| < 1, 

\^-v\< P, + V\ > P- (24) 
From the choice of ip and formulas above it follows, that 



Q<Z2-t + zi + + (21 + Z2Y < At. (25) 
While TT - 6», < 6'2 < TT, 

< ^ < e = tg9c. (26) 

Besides that, the inequalities < 92 < t — ti, < 9i < t holds authomatically 
because ti + T2 = t — zi < t as T2 > 0, and zi >= 1. Due to (PS)) . 

> zi > 
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Inequality (pS)) with fixed zi and Z2 gives restrictions for the variable ^ 

e„^^n = (t ^ 2zi)(t - 2zi) > > _ 2zi)(t - 2zi) + 



Cross-section of area of integration by the plane parallel to the plane zi = Z2 = 
^ = is presented on the Fig. 1 in the case of ^ < po, ez2 > ^ + po • It is seen 
that the straight line i] = €Z2 in dependence on Z2,^,e crosses shaded area or 
not. It determines the integral (IT9|) evaluation. 



5 Evaluation and approximate formulas 

Expression for l2{t) splits to three contributions hi , I22 , ha] h2,h3 respond 
the case ez2 < ^max + Pa- For the integral h2 Po < Cmax, while for the ha an 
opposite equality £,max < Po holds. 

One can show that in conditions of eOmaxPoln— << 1 and | >> where 
Cmaa; = max^ cr(cos7) , integrals h2 and ha give contributions small compared 
to one of one-fold scattering. In particuolar even rough estimation gives: 

h2 < , (.^yj 

ha < i^^^!4^(2 + ^ + l^^)- (30) 

The integral I12 in the same conditions is comparable with the contribution 
of one-fold scattering. Limits of integration in it are found from the condition 
£Z2 > ^rnax + Po- This couditiou determines the part of the square (j27p . points 
of which satisfy inequality 

Denoting by the simbol Z^At the integration domain of the plane Z1Z2, one has 
hi = limAt^o In^^ dz,dz2 C dp *y 

Ea{cos9\ ,z\)a{cos^ ,Z2)2p^ri V'^^J 
{(zt-Z2)'^+e){t-y-l-y/{7-l-7-2r+e){P^-((.-v)'')(p''-{(.+ri)^)' 



where E is defined above (l20l) , cos9i is given by (l22t . If one suppose that in a 
range of small cos^ (scattering angles close to ^) the cross-section a changes not 
very quickly, the expression <t{cos^) may be approximated by a{cos{'K — 9i)), 
or, as in 5 , by a{cos{'K — 9 — 9^/2)). 

In such conditions the integration with respect to 7y, p, ^ can be performed 
in explicit form. After transition to limit of Ai — >■ 0, cos9i = . and 
the integral under consideration takes the form 

J 022/,. Ea(-eos6>i,Z2)o-(cos6>i,zi) 

hi = 27r Po / dzidz2 ttt ^ , (33) 

J Do zf(t~Zl~Z2) 
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where Dq is defined by following inequalities 




(£^2 - pI) 
2(t~2z2)' 



(34) 



In variables ^i, Z2 the exponential factor is: 



E = exp[-2 




(Jt{T)dT 



t 



t-2zi 




at{T)dT 



(35) 



The expression p3p is the main result of the article: it represents, in the as- 
sumptions made here, the two-fold contribution into a response of the LIDAR 
probing. 

6 Conclusion 

The expansion we use converges because the number of particles conserwation 
law in elastic collisions. Estimation of of the convergence speed we delivered 
for a homogeneous layer thickness d. In this case the expansion by scattering 
multiplicity is the power series by optic thickness ad. 

Consideration of the three-fold contribution into a response, obtained by 
analogeous scheme in the case of Umaxd < 1 allows to establish conditions in 
which one can restrict himself by two-terms of the expansion. In this case, in 
principle, the inverse problem of the function a reconstruction may be posed. 

The polarisation account may be held by the scheme we use and do not 
introduce principal complications. 

A diffusion of the mean phase and interferention also may be included in 
this apparatus and do not contradict simplifications procedure. 

The theory development include more detailed modeling of the scattering 
processes and geometry conditions. 

By the way the condition amaxd < 1 is not necessary in a case of non- 
spherical indicatrix which is met in real LIDAR problems. The case of strongly 
nonspherical indicatrix, perhaps, can be considered on a way of unification of 
our results with ones of f6|. 

In this work we present relatively simple method for It means 
that we can control evaluation process. 

References 

[1] A.M. Kolchuzhkin, V.V.Uchaikin. Introduction into the Theory of Parti- 
cle Penetration through a Matter, Moscow, Atomizdat (1978) (in Russian). 
V.V.Uchaikin, A.A.Lagutin, The Stochastic Importance. Moscow, Energoat- 
omizdat (1993) (in Russian). V.V. Uchaikin, V.V. Ryzhov, The Stochastic 
Theory of High Energy Particle Transport, Novosibirsk, Nauka, Siberian 
Branch (1988) (in Russian). 



7 



[2] A. Buzdin, S. Lcblc. Boltzmann Equation Solution in LIDAR Problem The- 
ses Conference Lidar Zonding, Tomsk (1978) Theses Conference Lidar Zond- 
ing. Tomsk: 174-176, 1978. 

[3] A. Buzdin, S. Leble. Lidar Problem Solution in Double-Scattering Approxi- 
mation VINITI N2536-80dep p 126. 1980. 

[4] A. Buzdin, S. Leblc Boltzmann Equation Solution in LIDAR Problem Theses 
Conference Lidar Zonding,Tomsk (1978) Theses Conference Lidar Zonding. 
Tomsk: 174 176, 1978. 

[5] B. V. Kaul, I. V. Samochvalov. Equation of LIDAR probing of atmosphere 
with two-fold scattering account. Izv. VUZov, Fizika, 159, 109, 1975. I. V. 

Samokhvalov, "Double-scattering approximation of lidar equation for inho- 
mogeneous atmosphere," Opt. Lett. 4, 12-14 (1979). 

[6] M.V. Kazarnovskii, V. E. Pafomov, Solution of the Fokker-Plank equation 
for a layered madium. ZhETF, 74 , 846-848, 1978 



8 



